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1. Introduction 

Many different cosmological observations unanimously confirm the existence of a neutral, 
non-baryonic matter component in the universe. The dark matter density needed is 
five times that of ordinary matter and according to the standard model of cosmology 
constitutes about 20% of the energy content of the universe. However, current pieces of 
evidence stem from gravitational effects and the particle nature of the dark matter has 
yet to be established. 

Particle astrophysics is an active and exciting field but methods of indirect and 
direct detection of dark matter particles present in our surroundings suffer from 
uncertainties due to astrophysics and the unknown distribution of the dark matter. 
In high energy collisions at particle accelerators we can hope to create the dark matter 
particles and in a more controlled environment determine its properties. 

The CERN Large Hadron Collider (LHC) started up in November 2009. The plan 
is to let the proton beams collide at a center-of-mass energy of 7 TeV initially and if all 
goes well to ramp up the energy to 10 TeV, possibly via a series of intermediate steps 
with short periods of data taking. An integrated luminosity of at most a few hundred 
pb _1 is expected to be reached during this first run. 

Among the various new physics phenomena that may show up one of the most 
well studied hypotheses is weak scale super symmetry. Not only does low energy 
supersymmetry offer a solution to the hierarchy problem but this type of models typically 
also provide a perfectly stable TeV-scale neutralino - a particle which only interacts via 
the weak force. In general any such Weakly Interacting Massive Particle (WIMP) is a 
good dark matter candidate since it automatically leaves a Big Bang relic density of 
the correct order to explain the currently measured value, and the neutralino is not an 
exception - rather it is often used as a generic WIMP. 

At the LHC, heavy supersymmetric particles may be produced, the strongly 
interacting gluinos and squarks being the most abundant, and subsequently decay into 
lighter states via cascade decay chains. At every step of the decay chain ordinary 
standard model particles will be produced and give rise to reconstructed jets and leptons 
in the datajj] In general, isolation of a signal coming from new physics is far from trivial 
due to the large hadronic standard model background present. Hard cuts need to be 
imposed on the final state to discriminate signal from background. 

Many studies show that unlike the Higgs, superparticle states could be discovered 
early, i.e. at relatively low integrated luminosities. A stable WIMP would show up in 
the detector indirectly as missing transverse energy (MET), and search strategies often 
make use of combinations of cuts on MET and a number of jets and/or leptons with 
high transverse momentum (see e.g. [HISIE]). In principle the MET background from 
standard model physics arises from the production of neutrinos and should be quite 
moderate. 

Unfortunately, getting a handle on the MET measurement can be quite complicated 
\ In this context, "leptons" refers to electrons and muons only. 
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[3]. Just to mention some of the dominating contributions, dead and noisy channels in 
the calorimeters, beam-gas interactions, and pure mis-measurements in cracks, inactive 
material and escaping muons can fake large missing energy in the detector. Due to the 
above reasons it might not be feasible to make use of the MET channel when analyzing 
the earliest data. On the other hand, more locally measured objects like muons, 
electrons and jets can be identified without requiring as much global understanding 
of the detector as MET. This makes them less difficult to calibrate and more suitable 
for early measurements. 

Final states with many jets and leptons are rare in the standard model and as long 
as the jet and lepton multiplicities are large enough supersymmetry signal events should 
be visible above the background even without any requirements on MET, as shown in 
El [7] . Such events are hence ideal to search for even before any extensive calibration 
of the detectors have taken place. 

In this paper we study the prospects of early LHC detection of supersymmetric 
models interesting from a dark matter perspective. We impose no cuts on the MET but 
instead require high lepton and jet multiplicities within the events. Compared to the 
work in [5j|6j[7], where the constrained MSSM (or mSUGRA) was studied, we focus on 
the more phenomenological MSSM-7 supersymmetric framework, as well as an MSSM-8 
extension. We also impose the constraint that the models can provide the correct relic 
density of neutralinos to explain all dark matter in the universe. Another difference 
with our study is that we use a different set of tools (which have different systematics) . 

To make predictions using only Monte Carlo simulations before any data has been 
taken is difficult for various reasons that will be discussed below. Our aim at this stage 
is to estimate which parts of parameter space could be probed during early data taking 
without cutting on missing energy but otherwise using cuts that are fairly standard. 
Our predictions are valid for a center-of-mass energy of 10 TeV but we also discuss how 
our results change if a first data sample at 7 TeV would turn out to be useful. 

As knowledge of the detectors increases, other channels will become of high interest 
for supersymmetric dark matter studies. We wish to emphasize however that even 
after the detectors have been fully calibrated, channels without cuts on missing energy 
can provide valuable complementary information about the underlying physics model. 
This gives another reason for studying the prospects of detection of supersymmetry in 
channels not relying on missing energy. 

In Section 2 we give a short introduction to the dark matter problem. Section 3 
defines the supersymmetric model that we are using for our phenomenological study. 
The search strategy is described in Section 4 and the details of our analysis are given in 
Section 5. In Section 6 we present and discuss our results before concluding in Section 
7. 
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2. Dark Matter 

The evidence for dark matter in the universe is overwhelming (see e.g. [Hj) and suggests 
that the dark matter is made up by non-baryonic (elementary) particles. In particular, 
WMAP [9] has measured the dark matter density to be 

tt x h 2 = 0.1099 ±0.0062 (la), (1) 

where Q x is the dark matter density in units of the critical density and h is the Hubble 
parameter in units of 100 km s _1 Mpc -1 . In spite of the evidence, we still have not 
established the existence of dark matter other than through its gravitational interactions. 
There are however many hypothetical candidates for dark matter particles (see e.g. 
|10j ) and one of the most promising class of hypotheses is Weakly Interacting Massive 
Particles (WIMPs), that naturally (due to the weak scale interaction cross section) have 
a relic density of roughly the right magnitude as measured by e.g. WMAP, Eq. (CQ). A 
fascinating virtue of supersymmetric extensions of the standard model is that we get a 
WIMP dark matter candidate, the lightest neutralino, for free. Even if the neutralino 
naturally gives a relic density in the right ballpark, there is a large spread up and 
down depending on the details of the supersymmetric model. To be consistent with 
neutralinos making up all the dark matter of the Universe, we will therefore require that 
the lightest neutralino has a relic density within 2cr of the WMAP range, Eq. ([I]). This 
effectively excludes part of the supersymmetric parameter space, where the annihilation 
cross section is either too large or too small. 

Of course the neutralino does not need to make up all of the dark matter in the 
Universe, but it is much more appealing if we can explain the dark matter problem with 
just one new particle, especially one we get for free, as in supersymmetric models. One 
can also imagine that additional light states (that could arise in extensions of the MSSM 
models considered here) can affect the relic density of neutralinos via coannihilations. 
However, such light states will also affect the collider phenomenology studied here. To 
avoid these problems, we here stick to the simple scenario, where we assume that the 
neutralino makes up all of the dark matter, i.e. that the relic density satisfies Eq. ([Q). 

3. The Minimal Supersymmetric Standard Model 

The R-parity conserving MSSM contains 105 new parameters and it is not feasible to 
allow for all of these to vary freely. Hence, one usually makes simplifying assumptions 
based on various schemes for supersymmetry breaking. A common assumption is that 
of cMSSM (or mSUGRA), where gauge coupling and mass unification is assumed at the 
grand unified (GUT) scale and the parameters are then run down to the electroweak 
scale. 

We here choose to work with more phenomenological models, the MSSM-7 and 
MSSM-8, where we instead give the parameters at the electroweak scale directly. 
Compared to cMSSM, this gives us a larger parameter space as we essentially have 
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more freedom in the Higgs and neutralino/chargino sector. In the following subsections, 
we will briefly define our MSSM-7 and MSSM-8 models. 

3.1. Definition of the MSSM-7 

To reduce the 105 free parameters of the MSSM, a common assumption motivated 
by known low energy phenomenology is to set SUSY sources of CP violation to zero. 
Another is to assume the trilinear couplings and the soft sfermion mass matrices to be 
diagonal to avoid flavor changing neutral currents at tree level. 

For our MSSM model, we will first use the phenomenological MSSM-7 model, as 
defined in Ref. pTj . 

In this model, to reduce the number of parameters we make the approximation 
that only the third generation trilinear couplings A t and are non-zero and also 
introduce the common soft sfermion mass parameter m . In other words we take Ajj = 
diag(0, 0, A t ), A D = diag(0, 0, A b ), A E = and M Q = M v = M D = M L = M E = m l 
following the notation of [TT] . 

For simplicity, we also adopt the mSUGRA inspired GUT assumption for the 
gaugino mass parameters 

Q 

M 3 = — sin 2 6 W M 2 = -— cos 2 w Mi , (2) 
a 5 a 

where a and a s are the fine-structure and strong coupling constants, respectively. 
Approximately we have 1M\ ~ M 2 ~ O.3M3 and we choose to work with M 2 as our free 
parameter. 

Through the minimization of the potential for the Higgs scalar fields the mass 
parameters and the bilinear coupling in the Higgs sector can be reduced into two 
independent parameters: the already mentioned tan and the CP-odd Higgs boson 
mass itia- 

The seven free parameters left after these simplifying assumptions, fi, m^, tan/3, 
M 2 , mo, Ab and A t , defines the MSSM-7 reduction of the full MSSM parameter space 
pTj . Note that all model parameters of the MSSM-7 are to be given at the electroweak 
scale. This is different to mSUGRA where most parameters are instead specified at the 
GUT scale and run down to low energies using the renormalization group equations. 

3.2. Definition of the MSSM-8 

As a consequence of the assumptions made in the MSSM-7 above it follows that all 
sfermions, except for the top squarks and possibly the bottom squarks and tau sleptons, 
typically become nearly degenerate in mass. This feature is in fact more generic in the 
MSSM-7 than it is in for example mSUGRA, where the running of the masses down 
from the GUT scale can introduce extra divergences. 

While this mass degeneracy may be a convenient first assumption, it could also 
be overly restrictive when investigating properties that depend on the sfermion mass 
spectrum. A less degenerate spectrum opens up for new sparticle decay channels, and 
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hence for new signal characteristics. Since we in this work especially care about the 
production of leptons within SUSY decay chains, it makes sense to try to relax the 
relation between the squark and slepton masses somewhat. 

We therefore define the MSSM-8 as identical to the MSSM-7 except that we now 
split the common soft sfermion mass parameter mo into two: a common soft squark 
mass parameter rrig and a common soft slepton mass parameter mj. That is, we assume 
that Mq = Mj/ = Mj} = rrigl and = = mjl. In total we are left with eight 
free parameters: /x, rriA, tan/3, M 2 , rrig, mj, Af, and A t . 

4. SUSY search strategy 

In this work we set out to explore the prospects for detecting a signal from 
supersymmetry at the LHC, and the ATLAS detector [12] in particular, through 
search channels not relying on MET. Instead we focus on the signal of multijets plus 
multileptons that can arise in the sparticle decay chains characteristic of many SUSY 
models but is harder to produce within the standard model. 

With all the uncertainty that comes with a pre start-up simulation, the idea is to 
perform a quite conservative analysis. Specifically we keep our imposed sets of cuts 
simple and leave optimization until later, when the LHC has begun to collect some data 
reducing much of the uncertainty currently present within QCD phenomenology at these 
energies. 

For consistency we choose to work with a single computational package, 
MadGraph/MadEvent (MG/ME) [33] interfaced with Pythia [H] and the PGS [IS] 
detector simulation, for all signal and background processes. We make efforts to tune the 
parameters in our analysis chain so that our results agree with those presented by the 
ATLAS collaboration [16]. We also investigate sources of systematic errors by varying 
the settings around our default values. 

We analyze a couple of hundred models from wide scans over the SUSY 
parameter space, all required to pass numerous accelerator constraints and consistency 
requirements. An important restriction we impose on our models is that they have 
to provide the correct thermal relic abundance of neutralinos to account for all dark 
matter present in the universe. We calculate the neutralino relic density including 
coannihilations, threshold effects and resonances with DarkSUSY [T7]. Note that this 
approach is different from that of [5j [6j [7] where no constraint on the relic density is 
imposed but focus instead is put on specific regions of the mSUGRA parameter space. 

4-1- Search channels 

We choose a standard cut of four high energetic jets together with two isolated leptons, 
and furthermore divide the leptons into a same sign (SS) and an opposite sign (OS) 
search channel. In general both electrons and muons contribute to a dilepton signal but 
in case the identification of electrons within the early data turns out to be unreliable, 
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as argued in [BJ, we also perform an analysis in which only muons are counted. 

A standard procedure to gain control over the background is to require the 
reconstructed leptons to carry a transverse momentum, px, exceeding 20 GeV. It can 
however be argued that this requirement could be overly restrictive when it comes 
to muons, for which a cut on px > 10 GeV may be good enough [6j. Which value 
(or combination of different values) is optimal of course remains to be seen once the 
collection of real data has started, and for comparison we here simply choose to analyze 
both cases for the dimuon channel. 

Besides the dilepton search we also consider the channel of four jets plus three 
leptons. While [6j [7] found this channel to be slightly less interesting for the very early 
data, it can still be used as a cross-check. 

In order to investigate what impact the required jet multiplicity has, we also 
consider three-jet channels corresponding to the searches mentioned above. 

4-2. Signal properties 

R-parity conservation ensures that every sparticle decay always results in the creation 
of an odd number of new lighter sparticles. Two-body sparticle decays always produce 
one sparticle plus one standard model particle while three-body decays may give rise 
to one sparticle plus two ordinary particles. The new sparticle can then further decay 
into yet another lighter sparticle plus one or several ordinary particles, and so on. The 
decay chain finally ends with the perfectly stable LSP (which constitutes a dark matter 
candidate). The more steps in the chain, the more quarks and leptons can be produced 
and the clearer a signal we may get. 

At hadron colliders the dominant SUSY production channels typically are into 
color charged sparticles, i.e. into gluinos, g, and/or squarks, q. Consequently, the most 
important production processes at the LHC should be: pp — > gg, pp — > gq and pp — > qq 
(where q stands for any squark or anti-squark). 

The gluinos and squarks can then decay into lighter neutralinos, xl (i — 1)2,3,4), 
or charginos xt = 1)2) through processes like g — > qqx®., 9 ~^ Qq'xti Q ~~ QXa an d 
q q'xt- The neutralinos and charginos could further decay into lighter neutralinos, 
charginos or sleptons according to e.g. Xi ~~ Z^x°j, Xi ~^ ^xtji Xi ~^ W T ^xf, 
X° -> U, xt ~> W /±( * ) Xj > ) xt ~> hxf, xt ~^ Z^xf 01 X.t ~^ (where [ denotes any 
slepton or anti-slepton), with possible slepton decay back into neutralinos or charginos 
in accordance with I —¥ Ix® and I — > l'xf- Note that every one of the above mentioned 
decay processes can act as a source for jets and leptons reconstructed from the LHC 
data, either through quarks and charged leptons from direct emission or via the decay 
of the produced Higgs bosons, h, or gauge bosons, Z and W ± . 

The main factors determining the strength of a multijet plus multilepton signal 
are hence the gluino/squark production cross section as well as the structure of the 
spectrum of sparticles lighter than the gluinos/squarks. The production rate depends 
on the masses of the produced sparticles - the lighter the gluinos and/or the squarks are, 
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the higher the production cross section generally is. However, a sparticle spectrum with 
favourable masses and mass splittings could also help increase the signal. For example, 
the presence of light sleptons in the model should enhance the number of leptons emitted 
during the decay chain. 

With the above perspectives in mind we focus on models with light gluinos/squarks. 
To investigate the influence of the sparticle spectrum structure we analyze both MSSM-7 
and MSSM-8 models with various sparticle mass spectra. 

4-3. Background properties 

The standard model background can be severely constrained by demanding a large 
number of jets and leptons in the final state. The main background within the four 
jets plus two leptons channel comes from ti and Z production [7], for each of which we 
produce large inclusive samples with up to two explicit hard jets. Here, and throughout 
the rest of this paper, Z is used as shorthand notation for the full matrix element 
rf/zi*) ->. i+r. 

In the four jets plus three leptons channel the tiZ background may also be of 
importance [7] and we extensively simulate this process too. It should however be noted 
that due to practical restrictions on the maximum number of explicit jets feasible to 
include in the hard process within our analysis chain, our background simulation of 
the trilepton channel may not be as trustworthy as in the dilepton cases. (See Section 
16.21 ) This comes about since the dominating trilepton background seems to come from 
Z-events, for which an adequate number of explicit hard jets are hard to simulate. (See 
Section [6 . 1 . 1 j) . 

An important uncertainty in the determination of the background stems from the 
difficulty to correctly estimate the lepton fake rate, i.e. the amount of leptons not 
originating from electroweak processes like and Z decays. Fakes are typically 
produced in hadronization processes and do often propagate close to its parenting jet. 
This means that the fake rate is extremely sensitive to the imposed lepton isolation 
criteria (see Section 15. 3p . and due to uncertainties present in the modeling of QCD 
physics it cannot be determined to any greater precision by other means than from 
the upcoming data itself. In addition there are also true fakes arising from pure 
misidentification of non-leptons in the detector. 

Opposite sign lepton pairs are extensively produced in electroweak processes during 
ti and Z decays, and the OS dilepton search channel is not very much affected by the 
fake rate uncertainty. On the other hand same sign lepton pairs originating from ti 
and Z production must contain at least one fake. This implies that our simulated SS 
dilepton background is rather uncertain, and before the results in this channel can be 
trusted the lepton fake rate really has to be measured in the data. Trilepton events 
from ti and Z are of course also similarly affected by this fake rate uncertainty. 

We note however that our analysis chain seems capable of producing lepton fake 
rates of the same order as those estimated by the ATLAS collaboration in [T6], something 
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which adds some credibility to our SS dilepton background. For the pure SS dimuon 
channel, though, our fake rate seems to be roughly an order of magnitude below that of 

The main contributions to the background for the three-jet channels should also 
come from ti, Z and ttZ. 

5. Analysis 

We analyze the detection prospects for 100 MSSM-7 and 100 MSSM-8 models scattered 
over the (m squark ,m g [ uino )-p\d,ne\^ The models are extracted from extensive scans over 
the MSSM parameter space and all satisfy rriq < 1200 GeV and rrig < 1000 GeV, pass 
a large variety of accelerator constraints and provide a relic density of neutralino dark 
matter within la of the WMAP range, Eq. (OQ). For the MSSM-7, we have selected a set 
of typical representative models in various mass ranges, whereas for MSSM-8, we have 
focused on models with light sleptons. 

Within the MSSM-7, models with many different properties are chosen in order to 
search for correlations with the signal. Partly motivated by the discussion in [6], the 
mass splittings between the light squarks and the stops, sbottoms and staus, as well as 
between the neutralinos, are varied. Various values of tan (3 and fi are investigated. We 
also examine models where all neutralinos and charginos are lighter than the produced 
gluinos and squarks, as this gives many decay opportunities and hence possibly a larger 
signal. 

The chosen MSSM-8 models all include light sleptons, typically within the 100-200 
GeV mass range. In addition, we require all neutralinos and charginos to be lighter than 
the gluinos and squarks. These requirements open up for decays via the sleptons and 
hence a larger rate of production of final state leptons. 

For the signal we consider production of all possible two-body combinations of 
gluinos and ( ant i-) squarks. In total this adds up to 325 different hard process final 
states to simulate. We assume the production cross section into other states including 
neutralinos, charginos and sleptons to be sufficiently low to be safely neglected for our 
purposes. For each model we simulate 250 000 signal events. 

For the background, matched samples of ti + 0, 1, 2 jets, Z + 0, 1, 2 jets and tiZ are 
simulated and, in the cases of the ti and Z samples, normalized in agreement with NLO 
calculations [18]. In total we generate 12 • 10 6 ti events (~ 30 fb" 1 ), 26 • 10 6 Z events 
(~ 10 fb" 1 ), and 6 • 10 6 tiZ events (~ 10 5 fb^ 1 ). 

The signal and background samples are passed through a simplified model of the 
ATLAS detector in which jets and leptons are reconstructed, whereafter we impose cuts 
appropriate for our analysis. Main focus is put on the four-jet plus dilepton search 
channels but also trilepton as well as three-jet channels are examined. 

§ Here m squar k sets out to represent a generic squark mass scale, for which we take the up-squark 
mass m-a L (which is also typically close to the model input parameter m,q) as a good representative. 
Subsequent plots of the (TO sgua rfc, Wg/wino )-plane hence actually show the physical (m„ l , mg)-plane. 
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For each model we estimate the minimum integrated luminosity required for 
detection at the LHC by demanding 

S > max{5^B, 5}, (3) 

where S and B are the post-cuts number of expected signal and background events, 
respectively That is, we require a 5a significance but also at least five signal events. 

5.1. Tools 

The MSSM particle spectrum and relic density calculations are carried out with 
DarkSUSY-5 . . 4 [T7], a comprehensive numerical package for neutralino dark matter 
calculations. DarkSUSY also imposes a large variety of accelerator constraints on the 
models, including bounds on sparticle and higgs boson masses. The decay widths of the 
sparticles are calculated with SDECAY-1 . la [T9] . 

For the event generation of the standard model background as well as the 
supersymmetry signal we use MadGraph/MadEvent-4.4.21 (MG/ME) [13]. The MG/ME 
package includes a number of interfaced programs: MadGraphll identifies the relevant 
tree-level diagrams for the hard process and finds analytical expressions for the 
production cross section, while MadEvent constitutes the parton level Monte Carlo event 
generator. The events are then fed into Pythia-PGS-2 . . 27 via MG/ME's interface that 
includes the possibility of jet matching for generation of inclusive samples. Except 
for r-decays, which is handled by Tauola, Pythia takes care of the parton showering, 
cascade decays and hadronization. In the last step, the events are passed through PGS, 
the Pretty Good Simulation code which sets out to be a fast, but still fairly accurate, 
detector simulation. 

5.2. Settings 

We consider proton-proton collisions at a 10 TeV center-of-mass energy and choose 
cteq611 [20] for the parton distribution function. 

We make use of MG/ME's option of variable renormalization and factorization scales, 
and for our standard model background processes with explicit extra jets in the final 
states our choice of matching procedure is the if^-jet MLM scheme [2~T] . 

In Pythia string fragmentation according to the Lund model is used. We include 
initial and final state radiation but no multiple interactions. The jet measure cutoff 
used in the matching scheme is set to QCUT=25 GeV. 

For PGS we use parameters tuned to mimic the ATLAS detector, as provided in 
MadGraph/MadEvent-4.4.21. There are 81x63 calorimeter cells spanning the rj x 
(pseudorapidity x azimuthal angle) plane with Ar) = 0.1 and full coverage in 0. The 
electromagnetic and hadronic calorimeter resolutions are taken to be 10%/\/E and 
80% / \/E, respectively. A magnetic field of 2 T is assumed and tracking is covered for 
\rj\ < 2.5. Electron and muon coverage are \r)\ < 3.0 and \rj\ < 2.4, respectively. Jets are 
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reconstructed using a cone algorithm with cluster finder cone size AR = 0.4{jj] 

By comparing the lepton efficiencies for a W /± -sample using our analysis chain with 
those presented by the ATLAS collaboration [16] we find that an imposed additional 
efficiency, on top of what comes out from our detector simulation, can lead to an even 
better agreement. For this reason, we explicitly apply a muon efficiency of 0.9 for all \r)\ 
and an electron efficiency of 0.85 when \rj\ < 1 and 0.65 otherwise. 

5.3. Cuts 

After an event has passed through the detector simulation we apply our own cuts, both 
in order to better trust the result of the detector response and to gain in signal over 
background. 

We only consider isolated leptons, meaning that they could not easily have 
originated in jets. In PGS, electrons are isolated by default but not muons. The 
requirements on electrons are that the transverse calorimeter energy in an (3 x 3) cell 
grid around the electron, excluding the cell with the electron, in total has to be less 
than 10% of the electron's Et and that the summed px of tracks within a AR = 0.4 
cone around the electron, again excluding the electron itself, must be less than 5 GeV. 
Also, the ratio of the calorimeter cell energy to the px has to be between 0.5 and 1.5. 
For muons we require the summed px, excluding the muon itself, within a AR = 0.4 
cone around the muon to be less than 10 GeV. 

In addition we require any lepton (electron or muon) to be separated from the 
nearest jet by at least AR = 0.3 and to have a minimum p? of 20 GeV. In the search 
channel where we consider only muons and not electrons, we also let the minimum px be 

10 GeV following [61 [7]. While lowering the px cut could potentially cause more muons 
from leptonically decaying b-jets to be isolated, there is still much room for tighter 
isolation criteria at low pt, using both tracks and calorimeter cells, without significant 
loss of efficiency. Since an accurate prediction using such criteria is difficult to make 
without a control sample from real data, we choose to leave them out in the present 
analysis. 

We ignore electrons within 1.37 < |?7| < 1.52 since the detector response is unreliable 
in that region [16] . 

In order to reduce the Z-background we cut away events where any pair of same 
flavour, opposite sign leptons has an invariant mass in the range 75 GeV< m(l + l~) < 105 
GeV. 

In the four-jet channels we require the leading jet to carry a transverse momentum 
Pt > 100 GeV while each of the three sub-leading jets must satisfy pt > 50 GeV, 
following the procedure of many supersymmetry searches. In the three-jet channels we 
drop the leading jet and only require three jets with > 50 GeV each. 

For the channels with only muons instead of leptons, for the eventuality that 
electron identification cannot be trusted, we count the electrons reconstructed by the 

11 Ai? = y/( m - m ) 2 + Wi - h) 2 
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Figure 1. Total gluino/squark production cross sections for representative sets of (a) 
MSSM-7 and (b) MSSM-8 models. The various symbols correspond to different ranges 
for the value of the cross section. 



detector simulation as jets. 

Since we are considering the prospects of early discovery we do not apply a 
transverse sphericity cut, the effect of which would drown in jet systematics at this 
point. (See Section EES ) 

We do not apply any cut on missing transverse energy. 

6. Results 

To get the reader oriented with the models under study we plot their positions in the 
( m squark, m giuinoYv^^^ lY1 Figure [TJ We divide the points into two plots, one for the 
MSSM-7 set of models and one for the MSSM-8 set. For each model the summed 
cross section for pair production of any combination of gluinos and/or (anti-)squarks 
is indicated. As expected, the gluino/squark production cross section shows up as 
strongly anti-correlated with the gluino/squark masses. The results for the MSSM-7 
and the MSSM-8 are very similar - the slepton masses do not affect the production of 
gluinos and squarks. 

In addition to our wider MSSM scans we choose three pairs of benchmark models 
for a closer investigation of their properties. Each pair consists of one MSSM-7 and 
one MSSM-8 model that differ only by the value of the slepton mass parameter mj. 
The various pairs are otherwise spread out in the (m squar k, m 5 / uino )-plane and in the 
MSSM parameter space in general. The benchmark model parameters are listed in 
Table [T] together with some selected sparticle masses, the relic density and the total 
gluino/squark production cross section. Note that model B7 gives a too high relic 
abundance of neutralinos, and is hence excluded on cosmological grounds. This one 
model is included in our analysis despite this fact, in the purpose of better understanding 
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327 


377 


m xt 
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m xt 


236 
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379 


a x h 2 


0.12 


0.10 


3.8 


0.11 


0.11 


0.10 


O-prod (Pb) 


35.0 


8.6 


34 


.7 



Table 1. Benchmark models with derived sparticle masses, dark matter relic density 
and total gluino/squark production cross section. Each pair contains one MSSM-7 
and one MSSM-8 model for which only the slcpton mass parameter m ; ~ differs. All 
quantities with mass dimension arc given in units of GcV. 

the effect that light sleptons have for the LHC signal. 

By performing cuts on the simulated data, the background and signal cross sections 
for each search channel are determined. The standard model background cross sections 
are summarized in Table [2j together with the statistical Poisson uncertainty from our 
simulations. In the same way, the signal cross sections for our benchmark models are 
presented in Table [3j 

Using Eq. (|3]) we can then predict, for each MSSM model in our scan, the minimum 
amount of integrated luminosity needed for detection of a signal against the standard 
model background in the ATLAS detector for proton-proton collisions at a center-of- 
mass energy of 10 TeV. On doing so we neglect the uncertainty in the background 
but propagate the statistical errors in the signal cross section. The resulting statistical 
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Process 


ft (fb) 


Z(fb) 


UZ (10- 3 fb) 


Total (fb) 


1 , Vl £1 Tl Tl 1 


4j+leptons of p T > 20GeV 


qq 


11.2 ±0.6 


0.3±0.2 


144.3 ± 0.9 




11.6 ±0.6 


OS 


132 ±2 


7.4 ±0.9 


281 ± 1 




140 ±2 


31 


0.16 ±0.07 


0.5 ±0.2 


26.1 ±0.4 




0.7 ±0.3 


V^/llcLLlllt-l 


4j±muons of pr > 20GeV 


qq 


0.22 ±0.08 


(< 0.1) 


37.3 ±0.5 




0.26 ±0.08 


OS 


36 ±1 
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40 ± 1 


3/x 
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( , n fi Tiripl 
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qq 
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(< 0.1) 


49.6 ±0.5 




3.2 ±0.3 


OS 


51 ± 1 


7.4 ±0.9 


199 ±1 




59 ±2 


3/x 


0.09 ±0.05 


(< 0.1) 


9.6 ±0.2 




0.10 ±0.2 


V^llCLllllCl 


3j±leptons of p T > 20GeV 


qq 


45 ± 1 


1.6 ±0.4 


475 ±2 




47 ±1 


OS 


670 ±5 


72±3 


812 ±2 




743 ±5 


31 


1.7±0.2 


2.3 ±0.5 


95.3 ±0.7 




4.1 ±0.6 


( , V) SI Tl Tl £>1 


3j±muons of pr > 20GeV 




1.6 ±0.2 


(< 0.1) 


129.3 ±0.9 




1.7±0.2 


OS 


173 ±2 


45 ±2 


351 ± 1 




219 ±3 


3/i 


(< 0.006) 


(< 0.1) 


18.6 ±0.3 


(18 


.6 ±0.3) ■ 10- 3 


Channel 


3j±muons of pt > lOGeV 




14.9 ±0.7 


(< 0.1) 


168 ±1 




15.1 ±0.7 


05 


253 ±3 


73 ±3 


546 ±2 




326 ±4 


3/i 


0.09 ±0.05 


0.1 ±0.1 


37.3 ±0.5 




0.2 ±0.1 



Table 2. Background cross sections for different search channels. The indicated errors 
are statistical. In the cases where our simulations have not generated a single event 
that pass the cuts, we quote an upper limit (set by assuming one event, given the 
generated luminosity) within parentheses. 

errors for the derived integrated luminosities are typically small (< 20%). The reaches 
in various search channels are presented in Figures 1714121 and the results are discussed 
in more detail in Sections I6.1H6.3I 

A summary of the resulting detection luminosities in the four-jet channels can 
also be found in Figure [21 where the integrated luminosity is plotted against the total 
gluino/squark production cross section. Overall, we find that the prospect for early 
detection is better for the MSSM-8 models than for the MSSM-7 ones, something which 
points to the importance of light sleptons in the sparticle decay chain. The correlation 
between production cross section and detection luminosity is found to be strongest in 
the OS channel, somewhat weaker in the SS channel and much less prominent in the 
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Model 


A7 (fb) 


A8 (fb) 


B7 (fb) 


B8 (fb) 


C7 (fb) 


C8 (fb) 


Channel 


4j+leptons o: 


: p T > 20GeV 


SS 
OS 
31 


30 ±2 
87 ±4 
4.6 ±0.8 


63 ±3 
200 ±5 
24 ±2 


21 ±0.9 
61 ± 1 
5.0 ±0.4 


148 ±2 
256 ±3 
40 ±1 


30 ±2 
107 ±4 
3.6 ±0.7 


46 ±3 
266 ±6 
14 ±1 


Channel 


i 


:j±muons of pr > 20GeV 


SS 
OS 
3/i 


8± 1 
36 ±2 
0.4 ±0.2 


19 ±2 
101 ±4 
5.2 ±0.9 


5.8 ±0.4 

27 ± 1 
0.9 ±0.2 


45 ± 1 
110 ±2 
8.5 ±0.5 


10 ±1 
53 ±3 
1.2 ±0.4 


14 ±1 
145 ±5 
3.9 ±0.7 


Channel 


i 


:j±muons of pr > lOGeV 


SS 
OS 
3/i 


20 ±2 
72 ±3 
4.2 ±0.8 


33 ±2 
160 ±5 
14 ±1 


8.9 ±0.6 
42 ±1 
1.5 ±0.2 


57 ±1 
152 ±2. 
15.5 ±0.7 


15 ±1 
108 ±4 
4.0 ±0.7 


23 ±2 
302 ±7 
11 ± 1 


Channel 


3j+leptons o: 


: p T > 20GeV 


SS 
OS 
3/x 


65 ±3 
187 ±5 
11 ± 1 


142 ±5 
457 ±8 
54 ±3 


31 ± 1 
91 ±2 
8.0 ±0.5 


229 ±3 
408 ±4 
67 ±2 


46 ±3 
161 ±5 
7± 1 


76 ±3 
394 ±7 
22 ±2 


Channel 


3j+muons of pr > 20GeV 


SS 
OS 
3ji 


19 ±2 
73 ±3 
2.2 ±0.6 


41 ±2 
225 ±6 
13 ± 1 


8.6 ±0.5 
38 ±1 
1.4 ±0.2 


68 ±2 
169 ±2 
14.8 ±0.7 


15 ±1 
81 ±3 
1.5 ±0.5 


21 ±2 
213 ±5 
6.3 ±0.9 


Channel 


3j±muons of px > lOGeV 


SS 
OS 
3ji 


42 ±2 
152 ±5 
9± 1 


73 ±3 
347 ±7 
32 ±2 


12.7 ±0.7 

58 ± 1 
2.6 ±0.3 


86 ±2 
231 ±3 
27 ± 1 


23 ±2 
157 ±5 
5.8 ±0.9 


36 ±2 
434 ±8 
18 ±2 



Table 3. Signal cross sections for different search channels for the benchmark models 
defined in Table [TJ The indicated errors are statistical. 

trilepton channel. Regarding the other model parameters and sparticle mass splittings 
mentioned in Section [51 no clear-cut correlations with the signal magnitude could be 
established from our samples. 

The channel in which an individual model is first seen varies, although the SS 
dilepton channel typically looks like the most promising one. The possibility of using 
the real data to verify the background is however more solid in the OS dilepton channel, 
which may hence turn out to be of higher interest in upcoming analyses. Indeed the SS 
dilepton channel should not be much trusted at all before it has been measured in the 
data (see Section H~3l) . The results for the trilepton channel, which due to uncertainties 
in the background estimation should be taken only as rough guidance, are less interesting 
early on but should come into play as higher luminosities are reached. Often a model 
shows up in two or more channels, either at about the same time or successively at 
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relatively low luminosities still, offering the possibility to cross-correlate between them. 
6.1. The four-jet, dilepton channel 

6.1.1. Aj + 21 Figure [7] shows the resulting integrated luminosity required for detection 
in the four jets plus two leptons channel, to which both electrons and muons with 
Pt > 20 GeV contribute. The figure is divided into four parts showing the results for 
the SS and OS lepton cuts, for our MSSM-7 and MSSM-8 models respectively. Two main 
observations can be made: the MSSM-8 models have an advantage over the MSSM-7 
models in detectability, and the signal typically shows up earlier in the SS channel than 
in the OS channel. 

A majority of the examined MSSM-8 models with gluino masses below ~ 700 GeV 
turn out to be detectable within a few hundred pb" 1 of integrated luminosity. This is 
true for the signal in both the SS and the OS channel, although the latter typically shows 
up afterwords, and especially for high squark masses the advantage for early detection 
within the SS channel is apparent. 

Some MSSM-7 models with gluino masses up to ~ 600 GeV are detectable early 
on but generally only the region of lighter gluinos and squarks looks promising at low 
integrated luminosities. With a few exceptions, the SS channel seems to be the most 
promising one also here. 

The advantage of the SS channel may however be fragile. Our calculated 
background cross section, as given in Table EJ reveals that a luminosity of the order 
of 0.1 fb _1 must be collected before even a single SS background event is expected. In 
a real data analysis one would prefer to establish the background level by investigating 
the response to changes in the imposed cuts, for which a larger number of background 
events clearly is favorable. This is especially important for the SS background since it 
consists of lepton fake events (see Section H~3l) . The OS background, see Table EJ being 
an order of magnitude larger and consequently also easier to measure at low integrated 
luminosities, is on the other hand more robust. Hence the OS channel results may 
actually be the most important for the early searches. 

Another systematic uncertainty in the background determination stems from the 
fact that, for computational reasons, we include no more than two explicit jets in the 
hard process generation. This is principally a problem for the background coming from 
Z production, where ideally at least four explicit jets should be included. For the tt 
background two explicit jets should however give a good accuracy since e.g. t — > bW + 
and i —> bW~ decays can give rise to the remaining jets and isolated leptons. As we find 
the tt contribution to the background to clearly dominate the Z contribution within 
both the SS and OS channel (see Table [2]) we expect this kind of uncertainty to be 
moderate for the dilepton searches. 

6.1.2. 4j + 2fi In case electron identification does not become reliable until later stages, 
searches could still focus on muons. Figures [HH9] show the integrated luminosity required 
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Figure 2. Integrated luminosity required for detection versus total gluino/squark 
production cross section for the models in our (a),(c),(e) MSSM-7 and (b),(d),(f) 
MSSM-8 samples. Shown are the results for the (a),(b) SS dilepton, (c),(d) OS dilepton 
and (e),(f) trilepton search channels. 
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Figure 3. Integrated luminosity required for detection of the benchmark models of 
Table [T] in the dilepton (electrons+muons) and dimuon channels, at 10 TeV center-of- 
mass energy. The statistical errors (not shown in the figure) are less than 20%. 



for detection in the four jets plus two muons channels. Following the discussion in 
Section |4~T1 we show the results for two different cuts on the muon pt- 

Figure M shows the resulting detection luminosities for the dimuon search using the 
Pt > 20 GeV requirement. Overall the detection prospects are not as good as when 
electrons are included, although many MSSM-8 models with gluinos lighter than ~ 600 
GeV can still be found within the early data. Whether the SS or the OS channel is the 
most advantageous varies from model to model. The expected background in Table [2] 
indicates that several fb -1 must be collected before the first SS dimuon ti event shows 
up, which means that a real data verification of that background will hardly be feasible 
and the most solid early time results will probably rely on an OS dimuon analysis. 

In Figure [9] the corresponding results for the px > 10 GeV muon cut can be found. 
For basically every model, detection can happen sooner than for the higher muon p? 
requirement. For the OS channel the prospects are often even better than for the px > 20 
GeV dilepton search. Again the relative strength of the SS and OS signals varies, and 
after a couple of hundreds of pb" 1 many models with gluino masses up to ~ 400 GeV 
(MSSM-7) or ~ 600 GeV (MSSM-8) should be clearly detectable. 

6.1.3. Benchmark models In Figure [3] we show the resulting integrated luminosities 
needed for detection of our benchmark models, defined in Table (U in the four jets plus 
two leptons searches that we have considered. As the two models within each pair only 
differ in the slepton masses, the figure cleary illustrates the advantage of having sleptons 
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accessible in the sparticle decay chain. Note also that the relative significance of the 
various search channels varies between the models. 

6.2. The four-jet, trilepton channel 

As discussed in Section 16.11 the magnitude of especially the Z background comes with 
some uncertainty due to practical constraints on the number of explicit hard jets feasible 
to include in our simulation. Unlike what we find for the dilepton channels, Z production 
seems to dominate the background for the trilepton searches, see Table |2j This means 
that the systematic uncertainty may be rather severe within this channel. In addition, 
the sample size needed to find an accurate expected background is rather large, as can 
be seen from the statistical errors due to our simulation indicated in Table El 

Nevertheless, we have checked that our resulting background agrees fairly well with 
that found in [TJ, in particular in that the five-event level of Eq. ([3]) is what determines 
the detection up to fairly high integrated luminosities: ~ 1 fb" 1 for pr > 20 GeV leptons 
and ~ 10 fb" 1 for p? > 10 GeV muons in our case, see Table [2j 

We therefore, despite the uncertainty in the background determination, show plots 
of our resulting LHC reach within the > 20 GeV trilepton (electrons and muons) and 
the pt > 10 GeV trimuon channels in Figures fTOHTTl While detection of some MSSM-8 
models seems to be achievable before 200 pb _1 , the prospects are generally worse than 
for the dilepton channels, a conclusion which agrees with that of [6j [7J . The trilepton 
channels could be used later on to pin down the properties of the model underlying a 
detected signal. 

6.3. The three-jet channels 

Assuming the same processes to contribute to the background as in the four-jet searches, 
we show the results for the three jets plus two leptons channel in Figure [T2l A comparison 
with the corresponding four-jet results of Figure [7J indicates that the four-jet channel 
seems slightly more favorable, especially in the OS dilepton channel and for heavy 
gluinos and squarks. 

Basically, the same properties hold for the three-jet plus dimuon search, while the 
three-jet plus trilepton search suffers from large systematic errors within the simulated 
standard model background. Neither of those results are presented here. 

6.4- Systematic uncertainties 

The measure of significance of a signal in Eq. ([3]) assumes that the background 
uncertainty can be neglected. However, for the energies to be reached at the LHC 
a good knowledge of the standard model background cannot really be achieved before 
it can be measured in the data. As stated in the introduction, our aim with this 
paper is not to make exact predictions for very specific models but to get an idea of 
what supersymmetric dark matter models can be probed using early LHC data. There 
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■ pdf cteq6m 



minimum number of jets 



minimum number of jets 



Figure 4. Tests of jet systematics on the ti background in the dilepton channel. On 
the left we show the effect of a 20% uncertainty in the jet energy determination in the 
early data, and that the effect of a sphericity cut is small in comparison. On the right 
we show the effect of different choices for the renormalization and factorization scales 
as well as the effect of using a different pdf. The solid black line corresponds to our 
analysis. The scalefactor is an invoked rescaling factor common for all events. 



■ default cuts 

■ p_T-20% 
p_T+20% 

■ S_T>0.2 



- dynamical, scalefactor 1.0 

■ dynamical, scalefactor 0.5 
dynamical, scalefactor 1.5 

■ fixed at 2 mass 
fixed at gluino mass 



minimum number of jets 



minimum number of jets 



Figure 5. Tests of jet systematics on one of the benchmark models (B8) in the 
dilepton channel. On the left we show the effect of a 20% uncertainty in the jet 
energy determination in the early data. In the same plot, we show the effect of a 
sphericity cut on the signal. On the right we show the effect of different choices for 
the renormalization and factorization scales, where the solid black line corresponds to 
our analysis. The scalefactor is an invoked rescaling factor common for all events. 



are however some tests that we can perform on our simulated data to estimate how 
systematic uncertainties affect our results and to see how reasonable our cuts are. 

One difficult problem is that of what renormalization and factorization scales to 
use in the event generation. Additionally, for early data the systematic uncertainty in 
jet energy determination can be of the order 10-20%. 

Figure HI shows the effects of varying the jet pr cuts by ±20% (which closely imitates 
the effect of changing the jet energy) and the choice of renormalization/factorization 
scale for a ti sample in the dilepton channeljjj] Figure [5] shows the effects of systematics 

% A choice of a renormalization/factorization scale scale fixed at the Z mass is clearly questionable for 
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Figure 6. Integrated luminosity required for detection of the MSSM-8 benchmark 
models of Table [T] in the OS dilcpton (electrons+muons) and OS dimuon channels, at 
7 TeV as well as 10 TeV center-of-mass energy. 

for our benchmark model B8 correspondingly. (The systematics for all our benchmark 
models look similar). These systematic errors are not included in our analysis but should 
not alter the prediction for the detection luminosity by more than a factor of 2. 

We also show that the effect on the background of cutting away events with 
transverse sphericity less than 0.2 drowns in the uncertainty in jet energy determination 
within the early data, which is why we do not include this otherwise standard 
requirement in our cuts. We also note that the sphericity cut can reduce the signal 
at least as much as the background. 

We also note that while we here consistently assume a perfect knowledge of the 
true integrated luminosity, the experimentally measured value may not agree perfectly. 
This could also add some contribution to the systematics in the early data. 

6. 5. 7 TeV reach 

We investigate the effect of lower initial beam energy by running the simulation for our 
MSSM-8 benchmark models at 7 TeV center-of-mass energy. To get an estimate of the 
background we simulate ti events at 7 TeV to a luminosity of 2.3 fb -1 and focus on the 
dilepton channel for which we assume this to be the dominating background, and the 
OS channel i particular for which this gives a statistically representative sample. 

We see from Figure |6] that the integrated luminosity needed for detection at 7 TeV 

ti and, especially, heavy sparticle production, but is anyway included here for the sake of illustration. 
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is about 3-4 times higher than at 10 TeV. This is roughly what one would naively 
expect from how the production cross sections scale down, together with the effect of 
the decreased energy available for production of initial state QCD radiation. 

7. Summary and conclusions 

We have investigated early discovery signals in the ATLAS detector at the LHC without 
using a cut on missing transverse energy. Compared with the earlier studies of [SI [7J , 
we have instead of focusing on the constrained MSSM (or mSUGRA), studied the more 
phenomenological MSSM-7 and MSSM-8 models. We have also required that all our 
models have a neutralino that can be the bulk of the dark matter. In addition, we have 
used a completely different set of tools (with different systematics) than in previous 
studies. 

An important property to improve detectability within early LHC data, apart from 
lower gluino/squark masses, is the presence of light sleptons in the sparticle spectrum. 
In the MSSM-8 extension of the MSSM-7 model, we allow the slepton mass to be a free 
parameter and focus on models where sleptons are light. We conclude that a multi-jet 
plus dilepton search at 10 TeV, with no cut on MET, can probe supersymmetric models 
with gluino masses up to ~ 700 GeV for arbitrary squark masses (and even slightly 
heavier gluinos for lighter squarks) within the first 200-300 pb" 1 period of data taking. 
(See first and foremost Figured) The multi-jet plus trilepton channel is less important 
early on but can be used for cross-correlation at later stages. 

Our results seem to be fairly consistent with those of [SI El [7] regarding background 
and signal estimation, although a direct comparison is often not feasible due to 
differences in our analyses. We find our MSSM-7/MSSM-8 models to be a bit 
harder/easier to detect compared to their mSUGRA models under study, something 
which is consistent with the notion of the mSUGRA sfermion spectrum being a mixture 
of those of the MSSM-7 and MSSM-8. We note that even clearer signals may potentially 
be achieved by extending our MSSM-8 model to allow for separate mass parameters for 
the individual sleptons. 

It should be stressed that our imposed sets of cuts are intentionally very simple, 
and we expect that the LHC reach could be improved by performing more model 
specific searches. However, until the LHC has started collecting some data, providing 
a better understanding of the QCD background, large systematics remain and we find 
it reasonable not to go too much into detail. In fact, even though we have performed a 
conservative analysis in order not to overestimate the LHC reach, the systematic errors 
can still be quite important, and the results are to be seen as guidance. 

In summary, we have confirmed that supersymmetry searches, not relying on 
missing energy, using the early LHC data have a substantial reach within the MSSM 
parameter space, and found that models with light sleptons are particularly interesting 
to search for early on. 
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Figure 7. Integrated luminosity required for detection of (a),(c) MSSM-7 and (b),(d) 
MSSM-8 models in the four jets (with p T > 100, 50, 50, 50 GeV) plus two leptons (each 
withpr > 20 GeV) search. The results are shown for the (a),(b) same sign and (c),(d) 
opposite sign dilepton channels. 
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Figure 8. Integrated luminosity required for detection of (a),(c) MSSM-7 and (b),(d) 
MSSM-8 models in the four jets (with p T > 100, 50, 50, 50 GeV) plus two muons (each 
withpr > 20 GeV) search. The results are shown for the (a),(b) same sign and (c),(d) 
opposite sign dimuon channels. 
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Figure 9. Integrated luminosity required for detection of (a),(c) MSSM-7 and (b),(d) 
MSSM-8 models in the four jets (with p T > 100, 50, 50, 50 GeV) plus two muons (each 
wither > 10 GeV) search. The results are shown for the (a),(b) same sign and (c),(d) 
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Figure 10. Integrated luminosity required for detection of (a) MSSM-7 and (b) 
MSSM-8 models in the four jets (with p T > 100, 50, 50, 50 GeV) plus three leptons 
(each with pt > 20 GeV) search. 
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Figure 11. Integrated luminosity required for detection of (a) MSSM-7 and (b) 
MSSM-8 models in the four jets (with p T > 100, 50, 50, 50 GeV) plus three muons 
(each with pr > 10 GeV) search. 
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Figure 12. Integrated luminosity required for detection of (a),(c) MSSM-7 and (b),(d) 
MSSM-8 models in the three jets (with px > 50, 50, 50 GeV) plus two leptons (each 
with pr > 20 GeV) search. The results are shown for the (a),(b) same sign and (c),(d) 
opposite sign dilepton channels. 



